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Abstract 

A simple method to include the strong force in atom-antiatom scattering is pre- 
sented. It is based on the strong- force scattering length between the nucleon and 
antinucleon. Using this method elastic and annihilation cross sections are calculated 
for hydrogen-antihydrogen and helium-antihydrogen scattering. The results are 
compared to first-order perturbation theory using a pseudo potential. The pseudo- 
potential approach works fairly well for hydrogen-antihydrogen scattering, but fails 
for helium-antihydrogen scattering where strong- force effects are more prominent. 
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Antimatter in contact with matter causes annihilation - a violent process 
where mass is converted to energy. This is well known, even among people 
who never have studied physics antimatter annihilation seems to stimulate 
the imagination. In popular films and literature antimatter is used as e.g. an 
energy source, fuel for space ships, or a weapon. However, since the energy 
required to create antimatter outweighs the energy created by annihilation 
by many orders of magnitude, neither of these applications are possible in 
reality. More realistically, antimatter may have applications in medicine, such 
as antiproton treatment of cancer. 

The primary goal for antimatter research is not futuristic applications, but 
to uncover the symmetries (or possible asymmetries) of the fundamental laws 
of physics. One could expect that the matter and antimatter contents of the 
universe would be symmetric, but so far observations indicate that at least in 
the neighborhood of our galaxy there is only matter. Today the coexistence of 
matter and antimatter is ruled out on scales < 20 Mpc [1] . The reason for this 
matter-antimatter asymmetry could be a violation of the CPT theorem. This 
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theorem, and other matter-antimatter symmetries, could be tested through 
high-precision measurements on antihydrogen. Two experiments at CERN, 
ATHENA and ATRAP, have succeeded in creating antihydrogen atoms, but 
much work remains before CPT tests can be carried out [3] . 

Despite the interests in antimatter of both scientists and the general public, 
relatively little is known about the interaction between atoms made of anti- 
matter and ordinary atoms. So far no experimental data exist, and theoretical 
work has only just started. Atom-antiatom scattering has practical implica- 
tions for antihydrogen experiments, e.g. unwanted scattering of antiatoms on 
impurities or on the walls of the experiment, or scattering on ultracold atoms 
introduced on purpose to achieve sympathetic cooling of antimatter. Atom- 
antiatom scattering is also interesting in its own right, being a new fundamen- 
tal process involving both atomic and low-energy nuclear physics. 

A number of inelastic processes are possible in atom-antiatom scattering: 
nucleus- ant inucleus annihilation, electron-positron annihilation, rearrangement 
processes (such as formation of positronium), and even formation of metastable 
"molecules". The most important of these are nucleus-antinucleus annihila- 
tion and rearrangement. Rearrangement processes will also eventually lead 
to annihilation, but on a time scale much longer than the time required for 
the final-state fragments to separate. Hence, it can be regarded as a process 
distinct from the "direct" annihilation during the collision. Alternatively, re- 
arrangement can be regarded as a resonance in the annihilation channel. In 
this paper I will focus on the direct annihilation process, while rearrangement 
is discussed extensively elsewhere in these proceedings [4]. 

At large internuclear separations R the wave function of the atom-antiatom 
system with relative angular momentum / is given by 

$(R, n) -> sm(kiR - ln/2 + <S,)*im(^fl) ^(n^atom^i), (1) 

where an d a tom are the atomic wave functions of antihydrogen and the 
atom, ?*j denotes the coordinates of all leptons, and tiki is the relative momen- 
tum in the initial state. Both elastic and annihilation cross sections are given 
by the complex phase shift 5 t [2], 

a el = ^(2/+l) l-e 2t5 '\ (2) 
(7 - = ^(2/+l)(l- e - 4Im5i ). (3) 

At low temperatures (< 1 K) only the / = partial wave contributes, and the 
phase shift has the low-energy form 

lim tan5 = — ha, (4) 
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where a = a — if3 is the complex scattering length. Inserting (4) into (2) and 
(3) gives the low-energy form of the cross sections 



4tt (a 2 + (3 2 ) , 



(5) 



a a = —0. 
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We first note that the elastic cross section is constant at low energies, while 
the annihilation cross section diverges as k~ l . Thus, below a certain energy- 
annihilation dominates. Second, we note that <r el depends on both a and (3. 
This means that a non-zero imaginary part (3 of the scattering length, arising 
from the inclusion of the strong nuclear force, will not only give rise to anni- 
hilation, but will also significantly modify the elastic cross section. Only if the 
annihilation cross section is small, (3 <C \a\, this effect may be ignored. 

The simplest way to calculate the annihilation cross section is to use a pseudo 
potential 



This is a contact interaction with strength given by the constant A. This con- 
stant can be determined from experimental data on the nucleus-antinucleus 
system, or it can be calculated from model potentials. Using the pseudo poten- 
tial in first-order perturbation theory the annihilation cross section becomes 



For the proton-antiproton system A pp = 1.7 x 10~ 7 a.u.=6.8 x 10~ 37 eVm 3 , 
which was determined from the width 1130 eV of the Is state of protonium [5]. 
For the alpha particle-antiproton system no data on the Is state are available, 
instead the value A ap = 3.4 x 1(T 7 a.u.=1.4 x 10~ 36 eVm 3 was determined 
from a combination of low-energy annihilation data, the energy shift of the 
2p state, and the width of the 2p and 3d state [6]. Both values represent an 
average over different spin states. 

The pseudo-potential approach, being based on first-order perturbation the- 
ory, does not take into account the modification of the initial channel due 
to the strong force. Even if annihilation cannot be treated perturbatively one 
may still use the zero-range approximation of the strong force, since the atomic 
nucleus is very small on an atomic scale. The detailed shape of the strong- 
force potential will therefore not matter, only its strength as parameterized by 
a short-range strong- force phase shift 5 s f. Furthermore, atomic energies (eV) 
are very small compared to typical nuclear energies (MeV). The strong-force 
phase shift can therefore be parameterized by a strong-force scattering length 
a s f. This scattering length is complex, with the imaginary part representing 
annihilation. 




(7) 




(8) 
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The role of the strong-force scattering length can be understood by looking at 
the boundary condition of the wave function at short internuclear separations, 
R — > 0. In the absence of the strong force, the wave function has a short- 
distance form proportional to the regular Coulomb function for zero angular 
momentum F (kR). Even when the strong force is added the wave function 
is Coulombic for Rq < R <C lao, where Rq is just outside the range of the 
strong force. But now, since the Coulombic solution is not extended all the way 
to R — 0, also the irregular Coulomb function G (kR) is allowed. A general 
eigensolution to the zero-angular momentum Coulomb problem is an arbitrary 
linear combination of F (kR) and G (kR). For atom-antiatom scattering in 
the presence of the strong nuclear force the appropriate short-range form of 
the scattering function is given by a particular solution determined by the 
strong-force phase shift 5 s f 

= F (kR) + t a n5 sf G (kR) 
[ ^ l) F (kR )+t a n5 sf G (kR f [ l) - W 

Here </?(r*j) is the wave function of the leptons in the combined electric field of 
the overlapping nucleus and antinucleus. The wave vector k is related to ki as 
h 2 k 2 /{2n) = H 2 k 2 /(2pi) - E f , where E f is the leptonic energy corresponding 
to (f{r.j). Since Rq is very small on an atomic scale, the boundary condition 
of $ at Rq may be expressed in terms of the short-range expansions of the 
Coulomb functions [8], 

F (kR) -> C (v)kR, (10) 
G (kR) -> C^ir]) {1 - 2 V kR [\n(2rjkR) + h{rj) + 2 7 - 1]} . (11) 

Here the strength of the Coulomb interaction is given by rj = l/ib^k), where 
6 M = a m e / jj, is the mass-scaled Bohr radius. The normalization constant is 

C (r)) = \j2nr)/[l — exp(— 2rrr])], 7 ~ 0.5772 ... is the Euler constant, and h 
can be expressed in terms of the digamma function ^ as 

h(ri) = [V(-iri) + V(iri)]/2-hiri. (12) 

One then finds that [9] 

co t 5 sf - 2h( V ) = — ^, (13) 

where the strong-force scattering length is 

1 2tt . . 

= — — hm cot o s f. (14) 

a, f bf, fc-o 



The wave function $ may be expressed in terms of a regular and an irregular 
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solution as 

{It > l) [2h(r ] )-b,/a sf ]F (kR ) + C (r ] )y V G (kR ) ^ z) 
, 2h( V )F (kR ) + C%(ri)/riG (kRo) irrc 

+ [2/1(77) -bja sf ] F (kR ) + C (v) 2 /vGo(kRo) 1 ' 

(15) 

The regular and irregular solutions used are defined by their short-range forms 

1m Pi ) = FoikRMn) ~ ^(r 4 ), (16) 

R^o b^FotkRo) b^ 



lim $ irrcg (i2, r<) ! 



' |gMF (A; J Ro) + Go(A; J Ro) 



2r,h{v) F (kR) + G (kR) 



2fl 





<P(Vi 



+ 2 7 -l M^)- (17) 



These forms have the advantage that they are independent of Rq, provided 
that R is small enough. The normalization factors ensure that $ reg and <3> irreg 
are real functions, even if k is not. For large internuclear separations 

lim $ rcg (-R, Vi) = iV rcg 4 sm(kiR + 5 r c g )^H(n)0atom(n>i), (18) 
lim $ irrcg (#, n) = iV irrcg - sm(kiR + 5 irrC g)^H(^i)0atom(^i>i), (19) 

R—>oc Jx 

where both normalization factors and phase shifts are real. Inserting (18) and 
(19) into (15), and using (1), the complex phase shift Sq for atom-antiatom 
scattering including the strong nuclear force is found to be 

£ -^reg COS S TC g O s j/6^ t A^i rre g COS <5i r reg (20) 

iV reg sin 5 rcg - a s //^A^ irre g sin <5 irrcg 

Here the approximation kRo — > has been used. With this equation we arrive 
at a very simple way of calculating the strong-force effects beyond the pseudo 
potential in first-order perturbation theory. In addition to the usual regular 
solution (16), an additional irregular solution given by the boundary condition 
(17) has to be calculated. The asymptotic forms of the two solutions is then 
weighted according to (20) using the strong-force scattering length. 

The scattering length method has several advantages. It is very easy to im- 
plement numerically. It makes no assumptions about the detailed form of the 
strong-force potential, and is hence model independent. Alternatively one may 
directly use Eq. (9), in combination with Eq. (13), as a short-range boundary 
condition, and integrate towards R = oo. This was the procedure in Refs. 
[7,13,16]. The present formulation in terms of Eq. (20) has the advantage that 
it clearly brings out the relation between the atom-antiatom phase shift Sq 
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and the strong- force scattering length a s f, in terms of a minimal number of 
parameters. The parameters N Teg , N [rreg , 5 reg , and Si rreg are independent of the 
strong force. Once they have been determined from an atom-antiatom calcu- 
lation, it is a trivial task to calculate 5 for a given strong-force scattering 
length. 

Using this method elastic and annihilation cross sections have been calcu- 
lated for H-H and He-H. The strong-force scattering lengths used were a p ff = 
(0.844 - 0.698i) fm, and af f = (1.851 - 0.630i) fm, determined from the same 
data as the annihilation constants A pp and A ap above. Once again, these are 
averages over different spin states. The H-H potential was taken from Ref. 
[10] and the He-H potential from Ref. [11]. In figures 1 and 2 we compare 
the pseudo-potential and the scattering-length approaches for H-H and He-H 
respectively. 




Collision energy (atomic units) 



Fig. 1. Cross sections for H-H scattering; (solid) elastic cross section including 
the strong force, (dotted) elastic cross section without the strong force, (dashed) 
annihilation using the scattering length approach (dash-dotted) annihilation using 
pseudo potential. Only the s-wave part of the elastic cross section is shown. 

For H-H the pseudo-potential approach works quite well. As has been re- 
ported earlier, the low-energy annihilation cross section using the pseudo- 
potential approach is 0.14/^/qoq [12], while the scattering length approach 
gives 0.15/ v /qcIo [13] (e* is the collision energy). This is also in good agree- 
ment with a calculation by Armour et al. [14], using the optical potential of 
Kohno and Weise [15], who obtained / ^Jelal for the triplet state and 
0.15/y/elal for the singlet state. The elastic cross section changes from 771 al 
to 892 when the strong-force effects are included. The change obtained by 
Armour et al. [14] was from 788 al to 920 al for the triplet and 872 al for the 
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singlet. The scattering length is a = (8.4 — O.bli) a , so indeed (3 <C \a\ as 
required for the validity of the pseudo-potential approach. 

I 1 I 1 I 1 I I 




Collision energy (atomic units) 



Fig. 2. Cross sections for He-H scattering; (solid) elastic cross section including 
the strong force, (dotted) elastic cross section without the strong force, (dashed) 
annihilation using the scattering length approach (dash-dotted) annihilation using 
pseudo potential. 

For He-H scattering the situation is quite different. The strong force has 
a drastic influence on the elastic cross section. This is in part because the 
real part of the scattering length happens to have an untypically low value 
a = — 0.23ao (the typical scale of the scattering length is ~ 5 do). The elastic 
cross section without the strong force is therefore only 0.65 oo, but changes 
drastically to 1700 al when the strong force is included. This is reflected by the 
scattering length a = (2.2 — 11.4?) a , which has an imaginary part much larger 
than the real part. The annihilation cross section changes from 0.95/ \fela^ 
in the pseudo-potential approach to 2.62/ \feia1 in the scattering-length ap- 
proach. 

These results replace those in Ref. [16] which were based on the potential 
in Ref. [17]. This potential had been calculated at fewer R- values, giving in- 
terpolation problems in the cross section calculations. Rather surprisingly, 
the effect of this refinement of the potential on both elastic and annihilation 
cross sections was quite large. The scattering length reported in Ref. [16] was 
a = (—7.69 — 3.80i)a , giving an annihilation cross section 0.88/ \fela^. The 
increase of the calculated annihilation cross section with a factor three can 
almost entirely be attributed to the range R < 1.5ao of the potential. The 
difference between the old and the new potential is very small, at most about 
0.007 a.u. (at R = 0.8a ), showing the extreme sensitivity of the scattering 
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length to the details of the potential. Although this is a small correction, it is 
still 20 times larger than the maximum of the adiabatic correction [11], which 
on the contrary only has a small impact on the scattering length. 

The validity of the scattering length approach has been tested for the H-H 
system using a model potential with the form [6] 

K P t(R) = ~— (l + ^J 1 -) boP(R). (21) 
H \ A mj 

Here A is the nuclear mass number, b a complex parameter, \i the antiproton- 
nucleon reduced mass, m the mass of the nucleon, and p{R) the nuclear mass 
density, which was taken to have a Gaussian form with width tq. Three differ- 
ent widths were used, tq = 1.0, 1.5, 2.0 fm, and for each width the parameter 
b was adjusted until the potential reproduced the scattering length above. All 
three potentials thus obtained gave results identical to the scattering-length 
approach, to within three significant digits [7]. This clearly shows that the 
shape of the potential does not matter, and all its relevant properties can be 
summarized by the complex strong-force scattering length a s f. 

In conclusion I have developed a simple method to take the strong force into 
account in situations when the pseudo-potential method cannot be applied. 
This is the case in He-H scattering, where annihilation is fast, and the elas- 
tic cross section is drastically modified by the strong force. In the future this 
method may be used to investigate scaling properties of elastic and annihila- 
tion cross sections for scattering of antihydrogen on heavier atoms. 
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